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FINAL  REPORT 


1.  OBJECTIVES 

The  proposed  research  had  two  broad  objectives.  These  objectives,  as  stated  in  the  proposal, 
were: 

(i)  NETWORK  HYDRAULIC-GEOMETRY  AND  RUNOFF  PRODUCTION  To  seek 
analytical  generalizations  of  the  current  theory  of  channel  networks  which  include  hydraulic- 
geometric  variables,  e.g.,  channel  widths,  depths,  gradients,  velocities,  streamflows,  etc.  This 
will  require  theoredcal  studies  of  channel  networks  in  three-dimensional  space,  taking  into 
account  both  elevation  and  horizontal  spatial  dimmsims. 

(ii)  MESOSCALfe)  CLIMATE  To  undertake  analytical  and  empirical  studies  of 
intermittency  and  'statistical  self-similarity*  in  space-time  rainfall  intensity  at  natural  basin 
scales. 

2.  A  SUMMARY  OF  MAIN  FINDINGS  AND  RESULTS 

The  most  inqxirtant  theme  which  provided  the  focus  for  our  entire  research  effort  on  both 
of  the  above  objectives  is  an  understanding  of  the  scaling  invariance  property  in  spatial 
variability  of  channel  network  geometry,  rivor  runoff  and  rainfall.  We  will  first  explain  our 
main  findings  involving  scaling  invariance  on  both  of  these  objectives.  Then  as  an  application 
of  this  research  wt  will  briefly  describe  some  preliminaiy  results  cm  the  prediction  of  major 
floods  from  ungauged  basins.  A  solution  of  this  problem  is  important  to  the  Army,  and  is  the 
main  reason  for  the  support  ci  our  basic  research  by  ARO  since  1983 . 

OBJECTIVE  I 

Scaling  in  Channel  Network  Geometry 

Our  research  on  this  objective  was  directed  toward  a  mathematical  understanding  of  the 
empirical  main  channel  length-area  relationship  for  river  netwwks,  called  the  Hack’s  law, 
based  on  the  two  postulates  oi  the  randmn  model.  The  noost  significant  analytical  finding  was 
that  fmr  large  network  magnitude  m  the  avmge  main  channel  lengdi  L  varies  as  the  square 
root  of  m,  or  the  drainage  area  A  since  As2m-l,i.e., 

(1) 
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a 
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where  as  the  scale  factor,  denotes  the  mean  link  length.  It  has  been  rigorously  proved  for 
ccmstant  link  lengdis  [Waymire,  1989]  and  Ux  exponentially  distributed  link  lengths  [Gupta  et 
al.,  1990].  The  results  in  Gupta  et  al.  [1990]  suggested  that  Eq.  (1)  should  be  true  under  broad 
mathematical  ctmditions  on  link  length  probability  distributions  regardless  of  the  parametric 
f(»in  of  these  distributions.  Indeed,  for  the  class  of  distributions  having  finite  variance  this 
suggestion  was  verified  rigorously  in  the  article  by  Durrett,  Kesten,  and  Waymire  [1991].  In 
this  sense  Eq.  (1)  exhibits  universality.  This  result  also  solved  a  long  standing  problem  in 
theoretical  fluvial  geomcuphology. 

The  empirical  finding  of  Hack  however  shows  that  the  exponent  of  the  network  magnitude 
is  about  .57  rather  than  .5  as  predicted  above  by  the  random  nnodel.  An  explanation  of  this 
discrepancy  was  explored  by  us  in  different  directions.  For  example,  an  unpublished 
calculation  by  one  of  our  graduate  students,  Scott  Peckham,  based  on  the  empirical 
observation  that  the  average  meandering  frequency  and  amplitude  of  rivers  vary 
systematically  as  one  goes  from  the  drainage  divide  to  the  outlet,  suggests  an  exponent  larger 
than  .5  in  the  Hack's  law.  However,  a  theoretical  demcmstration  of  Hack's  law  in  charmel 
networks  requires  that  this  systematic  behavin’  in  the  river  geometry  be  translated  into  an 
appropriate  assumption  on  link  lengths.  We  think  that  this  assumption  is  tied  to  a  scaling 
invariance  property  in  link  lengths.  Since  scaling  invariance  is  me  the  major  focus  of  our 
research,  we  are  cmtinuing  to  explore  this  line  of  investigation  to  see  if  it  can  theoretically 
explain  the  observed  departure  from  .5  in  the  empirical  Hack's  law. 

In  a  recent  article  by  us  [Gupta  and  Waymire,  1989],  the  hypothesis  of  sinq)le  scaling 
behavior  in  link  heights  was  formulated  and  shown  to  give  good  predictions  of  the  empirical 
observations.  Ifowever,  the  same  predictions  also  suggested  a  need  for  generalization  of  the 
simple  scaling  hypothesis  in  link  heights  to  multiscaling.  The  notim  of  multiscaling  was 
introduced  by  Gupta  and  Waymire  [1990].  We  expanded  a  considerable  effort  in  obtaining  a 
closed  form  solution  oi  the  (non-linear)  prediction  equation  employed  in  testing  for  simple 
scaling  in  link  heights  [Waymire,  1991].  Even  though  this  equation  was  solved  numerically 
in  Gupta  and  Waymire  [1989],  the  numerical  approach  has  serious  limitations.  Consequendy, 
a  closed  form  or  semi-analytical  solution  is  required  for  carrying  out  tests  of  multiscaling  in 
link  heists.  Unfortunately,  the  article  by  Waymire  [1991]  on  this  problem  showed  that  the 
main  tools  (Laplace  transforms)  employed  in  solving  the  special  (linear)  case  of  tiiis  equation 
does  not  go  very  far  in  serving  the  general  nonlinear  equation.  Therefore,  this  problem  has 
turned  out  to  be  much  more  difficult  than  what  we  had  anticipated  earlier.  We  are  now 
exploring  other  avenues  to  test  for  multiscaling  in  link  heights. 
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Scaling  in  Annual  Peak  Flows 


River  basins  span  a  broad  range  of  spatial  scales  ranging  firom  about  10''  to  10‘  km^.  So 
from  the  view  point  of  predictions  from  ungauged  basins  it  is  necessary  to  understand  and 
model  the  spatial  variability  in  rainfall  and  river  flows  over  this  range  of  scales.  In  our  earlier 
work  [Cadavid,  1988]  it  was  found  that  the  instantaneous  annual  peak  flows  from  the 
Southeast  United  States  and  the  Appalachian  Mountains  do  not  obey  statistical  self-similarity 
or  simple  scaling.  In  fact,  it  was  first  shown  in  a  paper  by  Dawdy  [1961]  that  the  regional 
peak  flows  violate  the  assumption  of  index-flood  method,  which  is  tantamount  to  the 
demonstration  by  Cadavid  [1988]  that  they  do  not  obey  simple  scaling.  However,  theoretical 
research  on  regimial  flood  frequency  over  the  last  three  decades  has  not  been  able  to  abandon 
the  index-flood  method.  Therefore  it  has  not  had  much  impact  on  field  applicatitMis  in  so  far 
as  the  current  procedures  to  predicting  regimial  peak  flows  ctmdnue  to  be  dominated  by  pure 
empiricism.  Likewise,  arguments  have  been  made  against  self-similarity  in  spatial  rainfall  in 
the  recent  literature,  because  it  ccmtradicts  many  of  the  empirically  observed  features  in  data; 
see  e.g.,  Gupta  and  Waymire  [1992].  In  order  to  accomnoodate  the  empirical  scaling  features 
in  regional  peak  flows  and  spatial  rainfall,  we  are  exploring  the  new  theoretical  framework  of 
multiscaling  [Gupta  and  Waymire,  1990].  This  frameworit  is  currently  being  developed 
further  as  well  as  being  tested  by  us  peak  flow  observations  from  different  regitxis  of  the 
United  States. 

OBJECTIVE  n 

Scaling  in  Spatial  RainfoU 

A  physical  understanding  of  the  structure  of  annual  peak  flows  across  a  broad  range  of 
scales  is  intimately  tied  to  the  nature  of  spatial  variability  and  extremes  in  ^atial  rainfall. 
Theoretical  atten^ts  to  model  ^toe-time  rainfall  during  the  last  fifteen  years  have  evolved 
along  IXQ  separate  lines.  The  first  group  qjptoaches  is  based  on  an  assumed  hierarchy  of 
scales  in  spatial  rainfall  as  noted  by  numerous  empirical  interpretations  of  remotely  sensed 
observations.  The  second  group  approaches  rests  on  the  assunqitimi  (tf  self-sitnilarity,  or 
simple  scaling  in  the  probability  distributions  of  spatial  rainfall.  The  current  theoretical 
developments  ittvc^  a  common  modification  of  each  these  assunq>tions.  It  is  based  on 
the  notion  cH  spatial  random  truss  distribution  or  random  mul^ractal  measures  geiurated  by 
randcMn  cascades.  These  measures  are  highly  singular  and  can  ciqKure  both  the  scaling  as 
well  as  tile  extreme  variability  and  intermittency  observed  in  spatial  rainfall.  The  current 
matbenutical  foundations  of  tiie  theory  of  random  cascades  in  the  literature  have  mostly 
focused  on  the  ensemble  prtqierties  of  these  measures.  In  a  conqprehensive  article  we 
investigated  the  matiienutical  foundations  of  multifractal  measures  with  regard  to  their 

4 


spadat  sample  average  properties  [Holly  and  Waymiie,  1992].  In  this  paper  it  was  shown 
that  the  ensemble  properties  and  the  spatial  sample  average  properties  are  not  the  same  for 
random  cascades,  because  the  spatial  law  of  large  numbers  does  not  hold  due  to  strong 
conelatitxis.  Tests  of  these  results  were  carried  out  by  us  on  the  quux-time  rainfall  data  set 
from  GATE-I  [Gupta  and  Waymire,  1992].  In  particular  we  carried  out  the  space-time 
analysis  of  fractional  wetted  area  of  GATE-I  rainfall,  since  this  quantity  is  of  basic 
significance  in  understanding  the  production  of  runoff  in  space  as  well  as  in  a  variety  of  other 
hydrometeorologic  and  hydroclimatologic  studies.  These  preliminary  tests  show  that  the 
fractional  wetted  area  exhibits  the  spatial  scaling  as  predicted  by  the  theory.  This  analysis 
also  shows  evidence  of  the  breakdown  of  the  spatial  law  of  large  numbers  as  required  by  the 
theory.  These  findings  ate  quite  exiting  and  we  are  currently  exploring  testing  this  themy  on 
other  rainfall  data  sets. 

PRELIMINARY  RESULTS  ON  PREDICTION  OF  FLOODS  FROM  UNGAUGED  BASINS 

Suppose  that  peak  flows  obey  multiscaling  in  ^ace.  Let  the  scaling  exponent  of  the  n-th 
statistical  moment  be  denoted  by  n  s  2,*-.  Moreover,  suppose  that  link  slopes  in  a 

river  network  obey  statistical  self-similarity  with  a  scaling  exponent  0  [Gupta  and  Waymire, 
1989].  The  problem  is  to  predict  0  in  terms  of  the  scaling  exponents  describing  flows  and 

to  test  these  predictions  against  the  value  of  6  computed  directly  from  the  link  slope  data  in 
charmel  networks.  The  first  steps  towards  solving  this  problem  were  taken  by  Kapoor  [1990]. 
He  formulated  a  rxKioa  of  uniform  power  distribution  over  a  network  in  dynamic  equilibrium 
and  applied  this  ccvicept  to  predicting  die  exponoitO  as 

0-^2- Pi  (2) 

Using  flow  data  from  Brandywine  Creek  in  Pennsylvania  he  found  that  die  predicted  value  of 
0  lies  in  the  same  range  as  die  values  of  0  computed  by  Gupta  and  Waymire  [1989]  for  four 
river  basins;  the  value  of  0  for  die  Brandywine  Creek  was  not  available  to  him  at  the  time 
Kapoor  published  his  findings.  This  line  of  investigadon  also  provides  the  first  major  step 
towards  lining  network  geometry  with  channel  forming  peak  flows. 

This  line  of  investigation  will  play  an  important  rede  in  aolvinf  the  problem  of  pnedicrion 
from  unggyed  basins.  There  one  can  corxqnite  0  from  die  maps  and  can  dien  i^ly  Eq.  (2)  to 
predict  the  P's.  However,  diis  result  is  only  {neliminary  and  much  mme  foundational  work 
remains  to  be  done  on  solving  this  important  problem. 
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